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SUMMARY 


The  probability  that  within  a  future  large  second  sample  no 
failures  will  occur  before  the  expiration  of  a  safe  service  life 
estimated  from  a  small  first  sample  and  the  probability  that  the 
proportion  of  all  future  observations  failing  before  the  estimated 
safe  service  life  is  smaller  than  a  given  proportion,  are  the  two 
measures  of  safety  that  we  adopt  here. 

Assuming  the  logarithm  of  the  fatigue  life  is  normal  with 
known  variance,  we  derive  formulae  for  these  measures  of  safety. 
Setting  the  safe  life  as  some  fraction  of  the  mean  estimated  by 
the  first  sample,  we  then  compare  the  influence  of  other  parameters 
on  these  measures  of  safety. 

From  this  assumption  it  is  shown  that  one  has  virtually  as 
high  an  assurance  of  safety,  measured  by  the  first  criterion,  when 
using  only  the  minimum  of  the  first  sample, as  one  does  by  using  all 
the  observations  in  the  first  sample.  If  one  uses  the  standard 
second  criterion,  namely,  the  confidence  level  of  a  lower  tolerance 
bound,  as  a  measure  such  an  advantage  is  not  retained. 
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INTRODUCTION 


Suppose  we  have  a  structural  component  in  an  airplane  which  is 
subject  to ‘failure  from  fatigue,  perhaps  due  to  the  cyclic  loading  of 
the  ground-air-ground  cycle,  to  acoustic  loading,  or  to  any  other  well- 
defined  phenomenon.  We  ask  what  estimate  we  can  make  of  a  safe  service 
life,  or  more  specifically,  the  lower  quantities  of  the  distribution 
of  the  time  until  failure,  from  a  first  sample  such  as  the  simulated 
testing  of  a  small  number  of  specimens,  or  from  the  observation  of  the 
first  few  units  of  production. 

The  important  question  with  which  we  concern  ourselves  is,  what  is 
the  probability  that  the  weakest  component  within  a  fleet  of  airplanes 
to  be  produced,  or  coming  into  service,  will  fail  before  a  certain  time 
established  by  derating  the  life  estimated  from  the  first  sample.  The 
problem  is  to  determine  the  derating  procedure  so  there  i§  little  proba¬ 
bility  of  failure  within  the  future  fleet. 

We  let  X.  i-l,...,m  be  the  i^  failure  time  resulting  from  the 

l  ,m 

first  sample,  and  we  let  Y  i«l,.,.n  be  the  i^  failure  time  from 

l  ,n 

the  second  sample  of  components  in  service.  We  assume  that  the  initial 
sample  of  lives,  simulated  in  the  laboratory,  or  the  first  service  lives 
have  the  same  distribution  as  the  service  lives  later  on,  the  logarithm 
of  which  we  assume  is  normally  distributed.  We  label  the  common  distri¬ 


bution  of  life  with  F. 
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From  the  first  sample  taken,  we  have  the  data  X  ■  (X^  m) 

and  a  given  derating  function  d  which  we  use  to  obtain  the  derated  (safe) 
life  d(X).  The  probability  that  no  failures  will  occur  in  the  fleet  be¬ 
fore  the  safe  life  we  call  the  fleet  assurance,  labeled  a;  it  is  given  by 

(1.1)  '  a  -  P[Yl  n  >  d(&)]  -  |{l-F(d(x)])n  dG(x) 

X 

IV 

where  G  is  the  joint  distribution  of  the  sample  X. 

The  probability  that  a  proportion  of  at  least  3  of  all  future  ob¬ 
servations  will  not  fail  before  the  derated  life  we  call  the  confidence t 
labeled  y,  which  is 

(1.2)  y  -  P{l-F[d(X) ]  >  6}. 

Clearly  in  this  case  the  derated  life  is  just  a  lower  tolerance  bound  for 
the  population  described  by  the  distribution  F,  and  the  fleet  size  n 
does  not  enter  in. 

We  chose  the  function  d  to  be  some  fraction  of  an  estimate  of  the 
mean.  Our  problem  is  to  study  the  influence  of  various  parameters,  such  as 
sample  size  and  variance,  on  the  measures  of  safety,  while  keeping  in  mind 
the  economic  desirability  of  using  as  few  of  the  first  few  ordered  obser¬ 
vations  as  possible. 

2.  EVALUATION  OF  THE  SAFETY  FACTORS 
USING  FIRST  ORDERED  OBSERVATIONS 

We  assume  throughout  that  the  logarithm  of  the  fatigue  life  is  nor¬ 
mally  distributed  with  unknown  mean  y  and  known  standard  deviation  o, 
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and  the  derating  function  is  some  fraction  of  the  estimated  mean  life. 
Hence  from  (1.1)  we  have 

A 

(2.0.1)  a  *  P[Y  >  pv]  ,  for  some  0  <  p  <  1, 

1  |U 

2 

where  v  ■  expty  +  a  /2}  is  the  mean  (expected)  life  of  a  component  and 
v  is  the  estimate  of  v  made  from  X. 


A  current  procedure  is  to  take  a  small  sample  from  simulated  labora¬ 
tory  testing,  perhaps  m  -  2  or  m  *  3,  and  set 


(2.0.2) 


In  X,  . 
i,m 


As  we  shall  prove  later  on,  in  this  case,  we  have  for  the  two  safety 
indices  a,  Y 


(2.0.3) 


00 

a  m  f  +  ~J=  )  d0(x) 


(2.0.4)  y  *  $[ yin  U  -  Sg)3, 

where 

(2.0.4. 1)  -5-iS_E  +  £  *Ug)  -  6 


with  <P  the  cumulative  distribution 
defined  by 


function  of  the  standard  normal. 


(2.0.5) 


exp 


dt. 


Notice  that  neither  a  nor  y  depend  upon  y,  and  that  as 
m  grows  large  a  approaches  the  value  0n(£) ,  and  y  increases  to  the 

value  one  if  £  >  5  and  decreases  to  zero  if  £  <  4fi.  To  obtain  some 

p  p 

idea  of  the  magnitude  of  these  indices  we  take  n  m  200,  p  ■  t,»  0  m  ■  ^ 
(this  last  is  equivalent  to  taking  the  standard  deviation  of  the  common 

logarithm  of  life  to  be  equal  to  .2).  We  then  have  £  *  2.95,  and  from 
table  II  with  m  ■  3,  we  see  that  a  -  .64,  and  by  taking  4  *  2.326 

p 

(i.e.  6  «  .99)  we  can  compute  that  y  -  .86  from  tables  of  the  standard 

cumulative  normal  distribution. 

A  natural  question  to  ask  is:  If  we  increase  m  and  observe  only 
the  first  of  the  ordered  observations  (or  perhaps  the  first  few),  can  we 
increase  a  or  y?  If  the  answer  is  yes,  by  how  much? 

In  order  to  obtain  as  much  generality  as  possible,  we  consider  ob¬ 
serving  only  X.  ,...,X  where  r  may  be  1  up  to  m.  Thus  if  r  *  m, 
1  ,,m  r,m 

we  have  the  entire  sample.  If  r  ■  1  we  observe  only  the  life  of  the 
weakest  component  out  of  the  m  simulated  lives. 

Following  our  basic  assumption,  we  let 


THEOREM  1.  If  Z*  are  the  first  r  of  m  ordered  inde- 

•  J  'M  X  j  ffl 

p 

pendent  normal  observations  with  unknown  mean  y  but  known  variance  a  \ 

A 

the  maxim* m  likelihood  estimate  y  of  the  mean  is 
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(2.1) 


T 


r,  m 


where  x  ie  the  solution  of  [2.1.  3], 
rtm 


Proof.  The  first  r  of  m  ordered  normal  observations  have  the 
joint  density  on  -  ®  <  z^  <  • • •  <  z^  <  <*>• 


(2.1.1) 


ml 

(m-r) ! 


where  4>  is  the  density  of  (2.0.5). 


Then  letting  K  be  a  constant  independent  of  p,  and  L  be  the 
logarithm  of  the  Joint  density  (i.e.,  the  log-likelihood),  we  have 


(2.1.2) 


Letting 

tracting 


ip(x)  ■  4>(x)/$(x),  x  ■  (p-z^/o  and  by  adding  and  sub¬ 
in  the  summation  we  see  (2.1.2)  is  equivalent  with 


(2.1.3)  (m-r) 

We  call  the  unique  solution 
tonicity  of 


<Kx)  -  rx 


r 

I 

i-1 


z  -z , 
r  1 


x r  m  which  exists  by  the  decreasing  mono- 


•  • 


<*  it*.*  . 


If  the  graph  of  ( x )  »  as  given  in  figure  [1],  is  used,  an  approxi¬ 
mate  solution  to  (2.1.3)  can  be  readily  found.  The  tables  of  reference  [2], 
can  be  used  to  determine  the  root  more  accurately. 

Now  we  state  the 

A 

COROLLARY  1.  If  r  =  1  then  a  constant  depending  only 

upon  m  and  is  the  solution ,  call  it  x  ,  of  the  equation  in  t 

*( t )  =  J_ 

t  m-1 


and 


p 


+  o  •  x  . 
m 


We  present  a  table  of  values  of  t  in  table  I.  And  of  course,  we  can 

m 

obtain  the  result  of  (2.0.2): 


COROLLARY  2.  If  r  =  m,  then  x 

-  . —  m,m 


m 


—  2rz 

mrr 

i=l 


mo 


-Z,  J 


and 


. .  in 

p  =  Z  +  o  •  x  =  —  V  Z  . 

H  mtm  mtm  m  ** 

It  follows  from  (2.0.1)  that 

a  =  P[£nY1  >  £np  +  p  +  a^/2] 
i. ,  n 


(2.2.1)  a  -  PUn*1>m  >  *nP  +  zr>m  +  o  •  ir#m  +  a  12]. 


Letting  £  be  as  defined  in  (2. 0.4.1)  and  setting 


(2.2.2)  U^n  -  0  «.nY^n,  V^m  -  0  zifin  1  *  1»***»m 

we  have,  letting  A  (t)  be  the  left  hand  side  of  (2.1.3), 

r  ,m 


(2.2.3)  a  -  p|u,  >V  -C+A”1  |2(V  “v,  ) 

v  '  I  l,n  r,m  s  r,m  [**  r,m  i,m 


which  we  can  express  in  integral  form  as 


(2.3)  a-  f  /*n|5‘VAr!m  I^VVll  h<vl>‘“’Vr  ndvi 

J —oo Cyj  <  •  •  •  <y  <oo  J  * 


-oo<v  <  •  •  •  <v  <® 

1  r 


where  h(v^,...,vr)  is  the  joint  density  in  (2.1.1)  with  o  •  1. 

Unfortunately  this  integral  is  not  easily  simplified  in  the  general 
case.  But  to  obtain  some  idea  of  its  behavior,  we  give  the  explicit  ex¬ 
pressions  for  the  two  extreme  cases  r  ■  1  and  r  ■  m. 

For  r  ■  1  we  have  from  (2.2.3) 


a  -  P[U  n  >  V  -5+tJ 
1  ,n  1  ,m  m 


m 


oo 

J  (fiu-v- Tm)  [l-0(v)  ]m  1^(v)dv. 


By  using  the  fact  that  4)  is  even  we  obtain  a  function  of  £,  say, 


A1(C)  *  m  J  0n[C“Tm+x]0m  1(x)4»(x)dx. 


For  r  ■  m,  we  have  from  (2.2.3) 
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and  in  the  same  manner  as  above  we  have  a  function  of  S,  call  it 


00 

a2(0  ■  /  0n  +“~zU(x)dx. 
-®  v/m 1 


which  we  recognize  as  that  given  in  (2.0.3) 


The  comparison  of  the  values  of  and  for  n  =  200,  m  «  5, 

is  given  in  figure  2.  We  note  that  for  a  given  value  of  £  A^(£)  >  A 

for  £  <  3,  but  the  reverse  inequality  holds  for  £  >  3.  However,  there 
is  not  much  difference  considering  that  we  use  only  the  minimum  observation. 

A  tabulation  of  A2(£)  is  given  in  table  II  for  m  and  n.  As  a 
simplification  we  also  present  a  tabulation  not  of  A^(£)  for  various 
values  of  m  and  n  but  of 

00 

(2.4)  A3(p)  -  /  $n(p+x)d$m(x) 

—  00 

where  p  -  £-t  . 

m 

A  short  tabulation  of  A^Cp)  is  given  in  table  III.  This  accomplishes 
our  study  of  fleet  assurance. 

Now  we  turn  to  the  question  of  confidence  in  these  two  cases.  By 
definition  (1.2) 

A 

Y  -  P[F(pv)  <  1-0] 

where  F  is  here  the  log-normal  distribution  with  parameters  y  and  a2. 


Now 

F(PC)  -  fl/w  +  w  +  (JjH  -  vj  ■ 

Hence  by  letting  £  and  ;  be  as  defined  in  (2. 0.4.1) 

P 

A 

(2.5)  Y  -  P  (^  <  5  +  C^J 

-  P(Vr,m  +  ||  <Vr,*-Vl,m)l  '  5  +  'lJ 

This  can  also  be  expressed  as  an  integral  but  we  give  only  the  two  cases 
r  *  1,  r  «  m. 

For  r  «  1 

>  -  ‘ -  V  1  - 11  -*a  -  \  -  V1" 

(2.6)  Y  -  1  -0m(Tm  +  Cg  -  £) 

For  r  ■  m  we  obtain  the  result  stated  in  (2.0.4). 

(2.7)  Y  -  $(  /5«  -  Cg)] 

Again,  for  purposes  of  comparison,  we  Cake  the  two  functions 

Bi<5>  ■  1  -  -  e  +  ca> 

B2(£)  -  -  c6)]. 

which  we  present  in  figure  3  for  m  ■  0  -  .99. 

It  is  clear  that  I^U)  is  to  be  preferred,  it  being  the  greater  for 


almost  all  values  of  £. 
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3.  SAFETY  FACTORS  IN  THE  CASE  SOME  INITIAL 
OBSERVATIONS  ARE  LOST 

In  order  to  obtain  a  bit  more  generality  we  consider  the  eatimation 
problem  in  the  case  when  the  first  few  operations  may  not  have  been  ob¬ 
served.  This  could,  for  example,  be  the  case  when  the  phenomenon  under 
consideration  is  the  occurrence  of  a  fatigue  crack  of  a  specified  dimen¬ 
sion,  and  at  a  certain  time  regular  inspection  reveals  that  several  of 
these  have  already  occurred.  Hence  we  consider  observing  only 

X  . X  where  1  <  s  <  r  <  m.  Thus  if  s  ■  1,  r  *  m  we  have  the 

s,m  r,m 

entire  first  sample.  If  r  *  s  *  1  we  observe  only  the  life  of  the 
weakest  component  of  the  first  sample. 


THEOREM  2.  If  Z  , ...,Z  for  1  <s  <m  are  the  (r+1-8) 

—————  v  Sjffl  r,tf7 

ordered  observations  of  m  independent  normal  observations  with  unknown 

2 

mean  but  known  variance  a  t  the  maximum  likelihood  estimate  y  of  the 
mean  is 


y  =  Z  +  o'Tis^rjm) 
r,m  J  J 


where  T(str,m)  is  the  solution  of  equation  (2.1). 


Proof.  Proceeding  exactly  as  before  the  sth  through  the  rth 

ordered  observations  have  the  joint  density  on  °°<z  <  •••  <  z  <«> 

s  r 


m! 


(m-r) ! (s-1) ! 


0 


l  i-s 


Using  the  notation  of  Theorem  1  we  obtain 


<3-1)  0  )-  U-l)  *(“J— )  +  (m-r)  -  0. 

1"S 
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which  is  equivalent  with 


(3.2) 


(m-r)i^(T) 


(r-s+l)T 


We  call  the  solution  of  (3.2)  T(s,r,m),  which  exists  uniquely  since 

the  left  hand  side  is  a  decreasing  function  of  t.|| 

A  A 

Also,  note  that  T(l,r,m)  -  x 

r  ,m 


As  before,  the  graph  of  ^  as  given  in  figure  [1]  can  be  used  to 
obtain  an  approximate  solution  to  (3.2)  and  the  tables  of  reference  [2] 
can  be  utilized  to  obtain  the  root  more  accurately. 


The  problem  of  obtaining  a  simple  formula  for  the  assurance  or  confi¬ 
dence  in  this  general  case  seems  rather  difficult.  But  there  follows 
immediately  from  Theorem  2  the 

A 

COROLLARY  3.  If  a  =  r  then  T(rtrtm)  ia  a  constant  which  depends 
only  upon  r^rrij  call  it  t^  ^  which  is  the  solution  of  the  equation  in  t 

(m-r)ty(t)  -  (r-Vtyf-t)  -  t  =  0 


and 

A 

(3.3)  U  =  Z  +  o>t 

r,m  r,m 

Now  if  we  want  to  obtain  an  expression  for  the  assurance,  we  have 

A 

directly  from  (2.2.1)  by  replacing  t  by  t  and  using  the  notation 

r  ,m  r ,m 

of  (2.2.2) 


a 


P[U 


l,n 


-5+t 


+V  ] 
,m  r,m 
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Slnce  V  has  the  density  defined  by 
r,m 


r(r)^r  ^  r<fi(v) 


—oo  <  V  <  oo 


we  obtain 


(3.4)  a  ■  f  0n(C_tr  m+x)  Q*  1(_x)  0m  r(x)  <J>(x)dx. 

-'•OO 

A  table  of  this  integral  would  involve  three  parameters  so  it  would 
seem  more  practical  to  reserve  computation  for  those  cases  of  specific 
interest.  No  tabulation  is  made. 


For  the  confidence  in  this  case  we  have  from  equation  (2.5) 


Y  -  P[(Z  -Wt  -u)/a  <  £  +  C  ] 
r  ,m  r.m  1-p 


Y  -  P[V  <  5  -  t  +  5  ]. 

r  ,iu  r  ytn  i™p 


But  since  the  distribution  of  V  is  given  by 

r  ,m  6  J 


P  [ V  <  x]  =  1 

r,m 


r-1 

-  [1-0(x)  ]m  j 


m-r 

55  zL  ( j)  ^  (”x>  (x> 


we  can  obtain  values  of  the  confidence  y  by  using  first  tables  of  $ 
(reference  [2])  and  then  the  binomial  tables  (reference  [1])  to  calcu¬ 
late  it. 


We  call  attention  to  the  further  specialization  which  occurs  when 
the  inspection  reveals  only  the  strongest  item. 
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COROLLARY  4.  If  s  =  r  =  m  then  in  (3.3).  t 

_________  *  m  j  m 

was  defined,  in  Corollary  1  and  given  in  table  1. 


In  this  special  case 

00 

a  =  m  I  <f>n ( ^(x)<f>(x)dx 

^  fYl 


y  =  PW ; 


<  Z,+t.  +z,^  A 
m  s2-8 


-c  wTiere 
m 


i 


y  *Awvg. 
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4.  NUMERICAL  COMPARISONS 

We  make  the  e.i sumption  throughout  that  o  ■  .46.  This  is  equivalent 
to  taking  the  stai  dard  deviation  of  the  common  logarithm  of  the  life  to 
be  equal  to  .2. 

Let  us  see,  if  we  derate  by  fraction  p  *  and  we  ^ave  a  fleet 

size  of  n  ■  200,  whether  we  have  more  assurance  by  choosing  the  minimum 
of  five  observations  or  whether  we  have  more  assurance  by  using  all  three 
observations. 

Take  m  *  5,  r  *  1,  p  =  then  £  =  2.95,  and  for  n  ■  200 
we  find  from  table  I  that  A^(2.95)  =  .66,  which  is  about  twice  the  value 
of  A2<2.95)  for  m  =  3,  which  we  discussed  on  page  4. 

Suppose  we  want  the  assurance  to  be  .90  when  we  use  the  minimum  of 
five  observations.  We  ask  what  value  of  p  should  we  take?  From  table  III 
we  see  that  if  A3(p)  *  .90  then  p  *  2.6  since  p  =  £  -  t^,  £*3.66 

by  table  I.  Now  from  (2. 0.4.1)  we  have 


(4.1) 


p  =  exp  {-o£  -  o  /2} 


and  since  o  *  .46 


P  ■ 


-1.79 

e 


To  show  the  sensitivity  of  the  assurance  to  o,  let  us  suppose  that 
o  *  .2,  instead  of  .46.  From  equation  (4.1)  above  we  obtain 


P 


-.75 

e 


.47 


1 

2  ’ 
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The  tables  which  are  provided  are  of  course  very  sketchy  and  are 
included  only  to  give  an  idea  of  the  magnitude  of  the  functions  and 

A2  for  a  few  selected  values  of  n  and  m.  However,  the  programs  which 
we  used  to  compute  these  functions  are  available  by  writing  to  the  author, 
and  would-be  provided  upon  request. 
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0  12  3  4 

Figure  la 


Figure  lb 

Graphs  for  obtaining  the  solution  of  the  maximum  likelihood  estimate 
of  the  mean  when  only  the  weakest  component  is  observ  ' 
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Figure  2 


Graphical  comparison  of  the  assurance  that  no  failures  will  occur 
in  a  fleet  of  size  n  *  200,  using  only  the  minimum  value  out  of 
m  *  5  observations  for  A.  (£) ,  and  all  five  observations  for  A, 

4 


(5). 
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Figure  3 


Graphical  comparison  of  the  confidence  that  no  more  than  one  percent 
failures  will  occur  within  the  future  fleet  when  using  only  the  mini¬ 
mum  of  m  =  5  observations  for  B^S)  with  the  confidence  when  using 
all  five  observations  for 
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TABLE  I 


r 

m 

the  solution  of 

the  equation 

Mlm  i 

t  m-l 

m 

i 

m 

T 

m 

m 

2 

0.506 

16 

1.641 

3 

0.765 

18 

1.693 

4 

0.936 

20 

1.740 

5 

1.062 

22 

1.781 

6 

1.160 

24 

1.819 

7 

1.241 

26 

1.853 

8 

1.309 

28 

1.884 

9 

1.368 

30 

1.912 

10 

1.420 

35 

1.975 

11 

1.466 

40 

2.029 

12 

1.508 

45 

2.076 

13 

1.545 

50 

2.117 

14 

1.580 

‘  55 

2.153 

15 

1.611 

60 

2.187 

Interpolation  may  be  used  between  15  and  30  to  obtain  two-decimal- 
place  accuracy. 
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A  2(0 


TABLE  II 


[  +*p)d «x) 


n  «  200  n  -  200  n  =  200 

m  ■  3  m  =  5  m  =  10 


e 

a2(0 

1.00 

.005 

1.25 

.014 

1.50 

.034 

1.75 

.074 

2.00 

.142 

2.25 

.242 

2.50 

.369 

2.75 

.511 

3.00 

.649 

3.25 

.768 

3.50 

.859 

3.75 

.921 

4.00 

.959 

4.25 

.981 

4.50 

.991 

4.75 

.996 

5.00 

.999 

a2(0 

A,(C) 

.001 

.000 

.004 

.000 

.014 

.003 

.042 

.017 

.102 

.061 

.205 

.162 

.349 

.326 

.516 

.525 

.676 

.708 

.805 

.843 

.894 

.924 

.948 

.967 

.976 

.987 

.990 

.995 

.996 

.998 

.999 

.999 

.999 

1.000 

n  =  400  n  *  400  n  =•  400 

m  -  3  m  =  5  m  *  10 


€ 

A  2U) 

1.00 

.001 

1.25 

.005 

1.50 

.014 

1.75 

.035 

2.00 

.077 

2.25 

.147 

2.50 

.251 

2.75 

.382 

3.00 

.526 

3.25 

.666 

3.50 

.784 

3.75 

.872 

4.00 

.930 

4.25 

.965 

4.50 

.984 

4.75 

.993 

5.00 

.997 

a2(0 

a2(0 

.000 

.000 

.001 

.000 

.004 

.000 

.014 

.003 

.043 

.016 

.105 

.062 

.213 

.167 

.363 

.341 

.534 

.547 

.696 

.732 

.822 

.862 

.907 

.937 

.956 

.974 

.981 

.990 

.992 

.996 

.997 

.999 

.999 

1.000 

fir  if  f* 
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TABLE  III 

00 

A„ (p)  “  [  0n(p+x)d^n(x) 

^  ^-00 

n  -  400  n  *•  400  n  -  400 

m  =  10  m  *  20  m  =■  5 


p 

A^(p) 

A3(p) 

A3(p) 

0. 

0.024 

0.048 

0.012 

0.25 

0.049 

0.094 

0.025 

0.50 

0.092 

0.170 

0.048 

0.75 

0.160 

0.283 

0.086 

1.00 

0.260 

0.427 

0.145 

1.25 

0.388 

0.587 

0.229 

1.50 

0.533 

0.734 

0.337 

1.75 

0.675 

0.849 

0.463 

2.00 

0.797 

0.924 

0.595 

2.25 

0.886 

0.965 

0.718 

2.50 

0.942 

0.986 

0.820 

2.75 

0.974 

0.995 

0.895 

3.00 

0.989 

0.998 

0.944 

3.25 

0.996 

0.999 

0.973 

3.50 

0.998 

1.000 

0.988 

n  -  200  n  -  200  n  -  200 

111=  10  m  =  10  m  =  5 


P 

a3(p) 

a3(p) 

A3(p) 

0. 

0.048 

0.091 

0.024 

0.25 

0.089 

0.165 

0.047 

0.50 

0.156 

0.273 

0.083 

0.75 

0.252 

0.414 

0.141 

1.00 

0.376 

0.569 

0.222 

1.25 

0.517 

0.716 

0.327 

1.50 

0.658 

0.833 

0.447 

1.75 

0.781 

0.912 

0.580 

2.00 

0.873 

0.958 

0.701 

2.25 

0.933 

0.982 

0.806 

2.50 

0.968 

0.993 

0.883 

2.75 

0.986 

0.997 

0.937 

3.00 

0.994 

0.999 

0.967 

3.25 

0.998 

1.000 

0.985 

3.50 

0.999 

1.000 

0.992 

-23- 


REFERENCES 

1.  Tables  of  the  Binomial  Probability  Distribution.  National  Bureau 
of  Standards,  U.  S.  Government  (1950). 

2.  Tables  of  the  Normal  Probability  Integral,  the  Normal  Density  and 


Its  Normalized  Derivatives.  N.V.  Smirnov.  Macmillan  (1965). 


